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High pressure behaviour of liquid GeO^: a molecular dynamics study 
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High pressure behaviour of liquid GeC>2 is investigated by means of molecular dynamics simula- 
tions in the pressure range 0-20 GPa and at various temperatures. In agreement with the recent 
experiments (PRL, 92, 155506, 2004), Ge-O coordination increases under compression. For 1650 K, 
at ~11.5 GPa the structure becomes 6 coordinated with a discontinuous change in the density. The 
transition from the low density liquid (LDL) to high density liquid (HDL) is found to be reversible 
and the transition pressure increases with the temperature. At lower temperatures, the low density 
to high density transition is found to be continuous, but with the coexistence of two structures. 

PACS numbers: 64.70.Ja, 61.20.Ja, 02.70.Ns 
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I. INTRODUCTION 

Liquid SiC>2 is a major constituent of magma 3 , and its 
properties, in particular response to compression, has im- 
portant geophysical implications. Due to the similarities 
in the structures, physical and chemical properties of sil- 
ica and germania (Germanium dioxide, GeO^), the high 
pressure behaviour of the Ge02 has also been of consid- 
erable interest. In particular, several studies have shown 
that germania undergoes similar structural changes as 
silica, but at much lower pressures^, making it easier for 
experimental measurements. 

Recent experimental and theoretical studies have 
shown polymorphism in liquids i.e., these undergo liquid- 
liquid phase transitions^. For example, liquid carbon is 
predicted to undergo a liquid-liquid phase transition un- 
der pressure^ while such a transition has been observed 
in the case of liquid phosphorus^. A familiar network 
forming liquid, water, undergoes solid-state amorphous- 
amorphous transition^ and is suggested to have a liquid- 
liquid phase transition below its melting point 9 . Other 
liquids, such as Si, Ge, C, SiC>2 and Ge02, that have 
tctrahedral molecular bonding, also have possibilities for 
observing pressure induced polymorphism in the liquid 
phaseifi. As Ge02 is a prototype network forming liq- 
uid the high pressure studies of liquid Ge02 are quite 
important. 

Ge02 exists in two crystalline forms: a low density a- 
quartz (hexagonal) structure with tetrahcdral Ge-0 coor- 
dination and a high density rutile (tetragonal) structure 
with octahedral Ge-0 arrangement. The rutile structure 
is stable at ambient conditions, while a-quartz phase is 
stable above 1300 K^i. Silica also has similar phases: a- 
quartz and stishovite, latter being stable above 8 GPa. 
Both germania and silica form glass, with random net- 
work of corner sharing Ge04/Si04 tetrahedra. Upon 
compression both crystalline and amorphous Ge02 has 
been observed to undergo a structural phase transition 
from four coordinated Ge to six coordinated 4 . 

Recently, the high pressure behaviour of liquid 
Ge02 has been investigated experimentally^. With the 



help of in situ x-ray absorption measurements on liquid 
alkali germanate (Li20-4Ge02 composition) 12 , Ohtaka et 
oii inferred a sharp coordination change between 2.5-4 
GPa. In particular, the deduced Ge-0 bond distances 
exhibited a sharp increase in the pressure range men- 
tioned above and monotonically decreased thereafter. So 
as to understand the underlying structural and coordi- 
national changes in these conditions, it would be useful 
to carry out molecular dynamic simulations. A recent 
classical molecular dynamics simulation of liquid Ge02 2 
indicated a structural change in the between 4-8 GPa. 
In those simulations, carried out on a 576 atom cubic 
cell, the structures at different volumes were generated 
by scaling density, ranging from 3.0 g/cm 3 to 6.5 g/cm 3 
(corresponding to different pressures) and were subse- 
quently equilibrated at 3000 K and 1500 K. Aside from 
the fact that these (microcanonical ensemble) calcula- 
tions were done with a few hundred atoms, these did not 
show the sharp coordination change as seen in the exper- 
iment. Therefore, we have carried out a more detailed 
MD study, with a larger system, by compressing an ini- 
tial liquid structure in small steps (canonical ensemble) 
analogous to what is done in the experiments. 



II. METHOD 

For classical molecular dynamics (MD) simulations, 
the initial configuration was generated by melting Ge02 
crystal at 3000 K and then quenching the liquid to the 
required temperature. A macro cell consisting of 4374 
atoms (1458 Ge and 2916 O) was prepared by apply- 
ing symmetry operations on an a-Ge02 unit celU£. The 
system was then heated in steps of 500 K to 3000 K, al- 
lowing about 100 ps to equilibrate at each temperature. 
After the thermodynamic variables, viz, total energy, cell 
volume etc, had stabilized at 3000 K, temperature was 
lowered at a rate of 5 K/ps to 2000 K, 1650 K, 1500 K 
and 1200 K. 

Our simulations employed the inter-atomic pair poten- 
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TABLE I: Potential parameters 



Aij in kj/mol 


pij in A 


dj in kJ/(A e mol) 


Ge-O 7.9108 xlO b 


0.16315 


9000 


0-0 2.9259 xlO 5 


0.304404 


4985 


Effective charges 






gGc = 1.5, go=- 


0.75 





tials of Oeffner and Elliotlpia, i.e. Buckingham type, 
Vij(r) = — + A i:j exp (-r y /py) - % 

~ij r ij 

Their potential parameters, obtained by fitting to the 
energies obtained from the SCF calculations, are listed 
m Table HI These were derived, however, for crystalline 
GeC>2 and has been found to reproduce well the struc- 
tural and dynamical properties of both quartz and rutile 
phase. For our studies, parameters rescaled by Oeffner 
and Elliott to reproduce experimental vibrational density 
of states (VDOS) were chosen, while Gutierrez et al have 
used the unsealed ones. 

Simulations were performed in the constant Nose- 
Hoover NPT ensemble using the variable cell method^ 
as implemented in DL_POLY 16 . Verlet leapfrog integra- 
tion algorithm was employed with a time step of 2.0 fs. 
Periodic boundary conditions were used and the long- 
range electrostatic interactions were handled by Ewald 
summation methocUt 

Experimental melting temperature of germania is 
1378 K. But in MD simulations it may be different de- 
pending on the properties of the potentials used. To 
verify this, we heated crystalline GeC>2 in smaller steps 
(100 K) to about 2000 K allowing it to equilibrate com- 
pletely at each step. Around 1050 K a structural tran- 
sition to /3-quartz structure was observed with a den- 
sity drop of 2.9%. Another volume dicontinuity around 
1550 K resulting in a disordered phase, indicated melt- 
ing. Cooling from 3000 K, however, resulted in similar 
structures at all temperatures. 

After generating the liquid at the desired temperature, 
pressure was raised in steps of 1 GPa allowing the sys- 
tem about 100 ps at each pressure for complete equi- 
libration. Close to phase transitions equilibration was 
slow and took 250-300 ps. 



III. RESULTS 

Fig- HI shows that structures of liquid Ge02 at 1650 K 
undergoes substantial changes from 0.1 MPa to 15 GPa. 
P-V relationships for our simulations at various temper- 
atures are shown in Fig. [21 During compression volume 
exhibits discontinuity at 10 GPa for the simulation at 
1500 K, 11.5 GPa for 1650 K and at 17.5 GPa for 2000 
K. However, for 1200 K no discontinuous change is ob- 
served. 




FIG. 1: (Color online) Snapshot of (a) initial phase (0.1 MPa) 
and (b) high density phase (15 GPa). For clarity one Ge-0 
unit is shown separately for each phase. Big atoms are Ge 
and smaller ones are O. Only a small section from the centre 
of the super cell is shown. Both systems are at 1650 K. 
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FIG. 2: (Color online) PV diagram of liquid Ge0 2 at 2000 
K, 1500 K, 1650 K and 1200 K. P-V curve for the system 
at 1200 K varies continuously unlike other curves at higher 
temperatures. 



The observed volume collapse implies a structural 
phase transition, in agreement with the experimental ob- 
servations 1 . Our analysis shows that at 1200 K, cor- 
responding structural changes due to densification take 
place continuously. Upon release of pressure, systems re- 
turned to the original states at 3.5, 2.5 and 2 GPa for 
2000 K, 1650 K and 1500 K respectively. Interestingly, 
as explained later, the continuous changes observed by us 
at 1200 K are in agreement with those of Gutierrez et al 
in their simulations at 1500 K. Apart from the difference 
in density, system at 1500 K, 1650 K and 2000 K had 
similar structure in both low and high pressure phases. 
Hence, in the following sections we present detailed anal- 
yses of the changes at 1650 K and 1200 K only. 
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A. Simulations at 1650 K 

The low density phase of the liquid was composed 
of randomly connected Ge-0 tetrahedra as shown in 
Fig ^a)i& with a small fraction of 3 and 5 coordinated 
germanium atoms. Average Ge-0 bond length was found 
to be 1.81 A, which is close to the experimental value of 
1.82 Ai 

In Fig. |3| variation of Ge-0 bond length and average 
Ge-0 coordination with pressure are presented. Both are 
determined by finding the number of oxygen neighbours 
and their distances from each germanium atom in the 
simulation cell, within a sphere of radius equal to the first 
minimum in the radial distribution function, averaged 
over 1000 configurations. 
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FIG. 3: (Color online) Variation of average first neighbour 
Ge-O distance and average Ge-0 coordination number with 
pressure for simulations at 1650 K. 

In the experiments^, the deduced Ge-0 bond length 
was observed to reach a maximum (sal. 91 A at 3-4 
GPa) and gradually reduces thereafter. These observed 
changes were ascribed to a change in the Ge-0 coordina- 
tion from 4 to 6. Our studies confirm these conclusions, 
though the computed transition pressure is somewhat 
higher and the computed variation is more spread out 
in the pressure. However, in spite of a more gradual in- 
crease in both Ge-0 coordination and distance compared 
to the experiment, about 50% of change occurs between 
8-10 GPa, as seen from FigsEl&lU The differences may 
be due the fact that we have simulated pure Ge02, while 
the experiments are on alkali-gemanates. The fact that 
the melting temperature of this mixture is much less than 
that for the pure Ge02, implies that inter atomic inter- 
actions between Ge-O, and 0-0 atoms are different from 
that of pure germania. In addition, the difference be- 
tween the experiment and theory may in part be due to 
limitations of the pair-potentials used and also due to the 
non-hydrostatic stresses in the experiments. 



The experimental results of Ohtaka et al can not as- 
certain whether the transition state is composed of time 
and space averaged mixture of fourfold and sixfold states 
or includes intermediate states such as with the fivefold 
coordination. From our results shown in Fig. 0] it is clear 
that system goes through a five coordinated phase. 
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FIG. 4: (Color online) Variation of number of nearest oxy- 
gen neighbours for Ge atoms with pressure for simulations at 
1650 K. At intermediate pressures appreciable fraction of 5 
coordinated Ge-O units exist. 

X-ray structure factors for various pressures are pre- 
sented in Fig EI S{k) was determined by summing up 
Fourier transforms of pair correlation functions weighted 
by corresponding x-ray form factors. 

Note that at higher pressure S(k) shows more peaks 
suggesting that the corresponding phase may have some- 
what more intermediate range order than the parent 
phase. The simulated structure as shown in Fig PjTj) 
confirms this as the high density liquid has alternating 
layers of O and Ge similar to that in the rutile phase. 
However, within the layers atomic arrangement has only 
the short range order. 



B. Simulations at 1200 K 

As mentioned earlier, our studies at 1200 K gave sim- 
ilar results as obtained in the previous MD work carried 
out at 1500 K 2 . This may be partly due to the different 
potentials used and difference in methods adopted for 
simulation. They had carried out simulations on a mi- 
crocanonical system with density 3.65 g/cm 3 . Although 
initial and final structures are similar in terms of lack of 
order, we failed to observe any volume discontinuities for 
these simulations. Low pressure phase had a density of 
3.64 gm/cm 3 , average Ge-0 distance 1.79 A and average 
Ge-0 coordination as 4.02. When compressed, average 
Ge-0 distance increased to 1.88 A between 0-8 GPa and 
remained almost constant thereafter implying still evolv- 
ing structure. 
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TABLE II: Results 



Temperature (K) Transition pressure (GPa) Density (gm/cm 3 ) Ge-O coordination 

Compression Decompression at 0.1 MPa at 15 GPa at 0.1 MPa 15 GPa 

2000 VL5 3l5 3A5 EM 4U2 5.40 

1650 11.5 2.5 3.55 5.91 4.03 5.92 

1500 10.0 2.0 3.59 5.91 4.02 5.92 

1200 3.65 5.75 4.02 5.80 
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Ge-0 coordination number reached 5.74 by 20 GPa. 
Structure factors at higher pressures did not show any 
new peaks implying absence of any additional partial or- 
dering, unlike at higher temperatures. 

The first peak in Ge-Ge pair correlation function ini- 
tially located at 3.38 A was found to split into two at 
about 5 GPa, located at 2.72 A and 3.40 A, indicating 
that two kinds of Ge sites have formed. Comparing with 
results for the high density liquid phase, we see that Ge- 
Ge distance of 2.72 A corresponds to octahedral struc- 
ture, and hence it can be concluded that changes that 
take place abruptly at higher temperatures are sluggish 
here, and the initial tetrahedral structure has only par- 
tially transformed into octahedral structure at these pres- 
sures. 

IV. CONCLUSIONS 

From MD simulations of Ge02 liquid at various tem- 
peratures and in the pressure range 0-20 GPa, we observe 
that Ge02 liquid undergoes an abrupt first order phase 
transition from a low density tetrahedral structure to a 
high density octahedral structure at 1500 K, 1650 K and 
2000 K. At 1200 K this process is impeded due to re- 
duced kinetics and the system undergoes a continuous 
transition between these phases. Some of results from 
simulations arc summarised in Table HI] 



FIG. 5: (Color online) Structure factors at various pressures 
during compression for 1650 K. Inset shows partial structure 
factors S a /3(k) against k for GPa 
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